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Abstract. - We consider small systems of bosonic atoms rotating in a toroidal trap. Using 
the method of exact numerical diagonalization of the many-body Hamiltonian, we examine the 
transition from the Bose-Einstein condensed state to the Tonks-Girardeau state. The system 
supports persistent currents in a wide range between the two limits, even in the absence of 
Bose-Einstein condensation. 



Introduction. - The advances in experimental techniques in the field of cold atoms have 
made it possible to access experimentally strongly correlated states that had only been part 
of the mathematical apparatus of quantum theory up to now. A famous recent example is 
the Tonks-Girardeau limit of bosonic atoms [1] that has been realized in the two experimental 
studies of Refs. [2]. Confined in elongated, quasi one-dimensional traps - as it nowadays can 
be realized in many laboratories around the world - at low densities or strong interparticle 
interactions the boson gas exhibits some properties which resemble those of a fermionic system. 

Here, we consider interacting bosonic atoms which rotate in a toroidal trap. We examine 
the transition from the so-called Gross-Pitaevskii limit of weak interactions, where the gas 
forms a Bose-Einstein condensate, to the Tonks-Girardeau limit. In this transition, while 
Bose-Einstein condensation disappears, superfluidity persists in a much wider region. "Super- 
fluidity", throughout this paper, means the metastability of currents, or in other words, the 
existence of persistent currents [3-6] , which is defined more precisely below. 

The whole question of superfluidity is based upon two crucial issues: firstly, that there 
is an energy barrier between the current-carrying state and the current-free state. Secondly 
that no single-particle operator AV connects the two states, at least to some order. In 
this case the decay rate of the current due to some weak perturbation that dissipates angular 
momentum and energy, if necessary, is highly suppressed. In the present problem the criterion 
we use is that no single-particle operator connects the two states at least to first- order in this 
perturbation. As we will see, the two states are connected via /S.V to lower and lower order, 
as the system goes from the Gross-Pitaevskii to the Tonks-Girardeau regimes. As a result, the 
characteristic timescale of the decay decreases in the transition between these two limits. Still, 
a transition that is not allowed even to first-order, will take a very long timescale (compared 
to the lifetime of these gases) for sufficiently weak perturbations. This approach is equivalent 
to examining the dynamic structure factor, as in Refs. [3,4]. 
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It is known that Bose-Einstein condensation is neither a necessary, nor a sufficient condition 
for superfluidity [7]. The problem investigated here provides an example where: (i) for Bose- 
Einstein condensation to imply superfluidity, interactions are necessary, and (ii) superfluidity 
exists in a much wider region between the two limits, even in the absence of Bose-Einstein 
condensation. 

Assuming the usual contact interactions between the atoms, the many-body Hamiltonian 

is 

H = -^j V'^WVVW j V^(r)V^(r)V(r)V(r) dv. (1) 

Here M is the atom mass and C/q = 47r?i^a/M, where a is the scattering length for elastic 
atom-atom collisions. Since wc assume that the atoms arc confined in a toroidal trap, we 
express the bosonic annihilation operator ill in terms of the annihilation operators Cm of an 
atom with angular momentum mh, 'ip{r) = J2m CmS^"^^/W. Here V = 2'kRS is the volume 
of the torus with R its radius and S its cross section, with VS <C R. Measuring the energy 
in units of Eq = h'^/{2MR'^) we thus find [8,9] 

H = ^ l'^c\ci + I X] A4^rnCn h+l,m+n, (2) 

/ k,Lm,n 

where the dimensionless constant g equals AaR/ S. 

Transition between the Gross- Pitaevskii and the Tonks- Girardeau regimes. - In the 
Gross-Pitaevskii limit one has a Bose-Einstein condensate. In this case, the typical kinetic 
energy per atom is ?i^/(27\/i?^), and the typical interaction energy per atom is uVq. Here 
n = N/V is the atom density and N is the atom number. From these two energy scales, 
we see that in this limit g should not exceed ^fQP ~ 1/A''. In the opposite Tonks-Girardeau 
extreme, the typical kinetic energy per atom is on the order of h^kp/(2M) (just like the 
Fermi energy in a Fermi gas), with kp = o-k and a = N/ilirR) the density per unit length. 
The interaction energy is still on the order of nUo, and thus the characteristic value of u is 
of order a/S. Therefore, the coefficient aR/S is of order N in the Tonks-Girardeau limit, 
g > 9tg ^ N. The parameter 7 in the Lieb-Liniger model [8] is {-K/2)g/N. Then, for g N , 
7 is indeed of order unity. The above estimates are in agreement with the two higher curves 
of Fig. 2, as explained below. 

Gross-Pitaevskii limit and persistent currents. Although the results of the numerical 
diagonalization that we present below are much more general, it is instructive to start with a 
mean-field variational calculation that is valid in the Gross-Pitaevskii extreme. The existence 
of persistent currents can be seen easily in the following way [7,10, 11]. Let us start with the 
usual basis states of angular momentum mh, (j)mifi) = e^"^^ / V^ttR (normafized to A^), and 
consider a linear superposition of 4>o and 4>i, ipva,r{0) = Cq (/)o + Ci 01 = V 1 — ^ '/>o + e^'^\/l(j)i. 
Here (p is an arbitrary phase, and I is a parameter, being equal to the expectation value of the 
angular momentum per particle, I = L/N. For / = 1, V'var describes a vortex located at its 
center, while for Z = the vortex is at an infinite distance away from the torus. 

Because of the extra kinetic energy, the expectation value of the energy in the above state 
is higher for i = 1 than for Z = 0. The first decisive question is whether there exists an 
energy barrier for some intermediate value of I. Such a barrier is expected to be present for 
sufficiently strong repulsive interatomic interactions for the following reason: As I decreases, 
at some intermediate value between zero and unity, the vortex passes through the torus, which 
implies that a node forms in the atom density. This, however, costs interaction energy, since 
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for I — and I = 1 the density is homogeneous, which is the configuration with the lowest 
interaction energy. Therefore, for sufficiently strong interactions, such a barrier is expected 
to develop. 

Quantitatively, the expectation value of the Hamiltonian in the state tpva.r above is the 
Gross-Pitaevskii energy Eqp, which is (in units of Eq) 



with S — {N—l)g [12]. According to Eq. (Q, Eqp develops a local maximum at I = (1+5)/ {2S). 
In order for this to appear at a value of I smaller than unity, S > I, or g > 1/ (TV— 1). Therefore, 
this is the critical value of S (and of g) above which an energy barrier develops in the dispersion 
relation Eqp = Eqp {I). Equation © also coincides with the excitation spectrum calculated 
by Lieb [9] in the limit N oo and R ^ oo, with N/{2ttR) = a, finite, and for small 7, i.e., 
Eqp/N - Airj = 87r7(/ - P) in the units of energy of Refs. [8,9], h'^a^/2M. 

In order for our argument for the persistent currents to be complete, we need to consider the 
possibility of connecting the many-body states ^'o with L = and ^'at with L = N via some 
single-particle operator AV of the form AV = VoS{9). Only one Fourier component of the 
delta function connects the two states and provides the difference in the angular momentum 
between them, thus 



The above operator is precisely the one that enters in the calculation of the dynamic structure 
factor S{k,uj) for k = 2kF (or L = in our model). In Ref. [4] it has been argued that the 
vanishing of S{k = 2kF,tu = 0) implies that the decay rate of persistent currents is suppressed, 
in agreement with the present study, as we discuss below. 



In the Gross-Pitaevskii limit 1 5*^=0 1 is independent of 9j, and the matrix element vanishes; 
only high-order perturbation theory in AV then connects the states \'^l=o) and \'i'L=N)- 

Exact results from numerical diagonalization of the Hamiltonian. - To go beyond the 
mean- field approximation (which breaks down far away from the Gross-Pitaevskii limit of weak 
interactions), we have performed numerical diagonalization of the full many-body Hamiltonian 
H. We thus evaluate numerically the eigenstates of H, which are also eigenstates of the 
momentum operator L, and of the number operator N, for various values of g. 

Although "exact" (apart from the naturally necessary truncation of the Hilbert space), 
this method is quite restrictive in the values of A^, L, and g that can be considered. More 
specifically, the number of the basis states that need to be included increases as one approaches 
the Tonks- Girardeau regime. In the Gross-Pitaevskii limit for a non-rotating cloud the domi- 
nant contribution to the many-body state comes from the m — single-particle state. In the 
opposite extreme of the Tonks- Girardeau limit, even for a non-rotating cloud, the dominant 
contribution comes from the single-particle states with |m| < mmax ^ A^/2, since in this limit 
the many-body wavefunction has to build correlations of characteristic length R/N . Further- 
more, for a given interaction strength Uq and a given density n, h?' rn^^^/ (2AI R^) ^ tlUq, or 
m^a^x ~ ^9- Thus, more generally, mmax ~ max(A^/2, ^Ng). 



Eqp/N - 5/2 = (1 + 5)1 - 5f + 0{g), 



(3) 



N 




(4) 
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Fig. 1 - The points show the lowest eigenenergies of the Hamiltonian £ — £{L) for L = 0, . . . , 5 and 
N = 5. Also g = 0.1 (top), g — 1 (middle), and g = 10 (bottom). The curves show the energy 
as given by the mean-field approximation Egp [Eq. @] in the top graph and Etg [Eq. ®] in the 
bottom graph. 



In Fig. 1 we consider N = 5 atoms, and L = 0, . . . , 5, with rrimax — 20, which gives a 
change of the lowest eigenenergy that is less than 1% between rrimax — 19 and m„iax = 20, for 
the highest value of g = 10 that we used. For N = 5, gop ^ ^/N = 0.2, while gTG ^ N = 5. 
Figure 1 shows the lowest eigenenergy £ of the Hamiltonian as function of the total angular 
momentum L, for g = 0.1, g = 1, and g = 10. As argued above, for g = 0.1 the system is 
essentially in the Gross-Pitacvskii limit. Furthermore, there is no maximtim in £ = £{L), since 
the interaction is not strong enough (the mean-field approximation predicts that the critical 
value of g is 0.25 for 5 atoms). The dashed curve in the upper graph of Fig. 1 shows the result 
of the mean-field approximation. Clearly, it is possible to improve the mean-field energies 
by including states with higher \m\ in V'var- For g = 1, there is a maximum in £ = £{L) at 
L — 4. Finally, for g = 10, the system is closer to the Tonks-Girardeau limit, and there is 
still a maximum in £ — £{L). The dashed curve in the lower graph of Fig. 1 is the energy 
predicted in the Tonks-Girardeau limit [9] , 

Etg = - + NL, (6) 

plus corrections of order N. Such corrections are seen clearly in this graph. For L = 
the leading-order term in the energy is /12, which is the total energy of a Fermi gas in 
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Fig. 2 - Two higher curves: the function 1 — g^^\0), for N = 3 and = 6 atoms. Two fower curves: 
the matrix element |(*L=o|Ay|*i=jv>/Vb|^ for 7^ = 3 and N = 6 atoms. 

one dimension, NEp/S (with Ep the Fermi energy), for a value of the Fermi momentum 
kp = = N/{2R). The expression for Etq implies that the energy barrier in f = £{L) [i.e., 
the difference between the maximum in £{L) and £{L = N)] scales quadratically with A'' in 
the Tonks-Girardeau limit, in agreement with our numerical results. On the other hand, in 
the Gross-Pitaevskii limit, the maximum of the energy barrier scales linearly with N . This is 
the reason for the emergence of this barrier as one goes from the one limit to the other. 

An exact result. - The eigenenergies calculated above satisfy an exact formula, valid for 
any value of g, which is in agreement with the study of Felix Bloch [13]. More specifically, 
£{L') — £{L) = L' — L, which also implies that £{L = N) — £{L = 0) = N, or more generally 
£{1 = qN) - £{L = 0) = q^N, g = 0, 1, 2, . . . Geometrically £{L') - £{L) = L' - L implies 
that all the lines connecting the points with L and L' = N — L are parallel. 

To see this result, we recall that any state \L) is related to the state \L' = N — L) via some 
collective excitation of the system. The state \L) is a linear superposition of Fock states of the 
form K-m)^-", . . . , 0^«, . . . , m^™), with Y.m = rnN^n = L. Then, the state 

that consists of the Fock states | (-m + 1)^" ,...,(- O^i , l^o , ... , m^-^+i , (m + l)^-") 
with the same amplitudes as \L) is \L' = N—L), since it has A'' atoms, L' = X]m(~™+l)-^"i = 
N~L, and a total energy which is higher than that of \L) by ^^^J(— m + 1)^ — m^JiVm = L' — L 
purely because of the kinetic energy; the interaction energy is the same in \L) and since 
all the matrix elements of the interaction are the same, independent of m. Similar arguments 
apply in all the branches with pN < L,L' < {p + 1)N, with p= 1, 2, 3, . . . 

Persistent currents close to the Tonks-Girardeau limit. In this last part of our study we 
examine the lifetime of the current-carrying state with L = N close to the Tonks-Girardeau 
limit. Due to the Bose-Fermi mapping [1], for 7 — > oo the matrix element (5'L=o|AF|^'i=Ar) 
is identical to the one for non-interacting fermions and equal to Vq. The two lower ciirves 
in Fig. 2 show \{^ l=o\AV\^ l=n) /Vo\'^ , for = 3 and N = 6, which indeed tend to unity 
for large g. Therefore, for 7 ^ 00 there is a single-particle operator that connects the two 
states and destroys the persistent ciirrcnt for L = N . In addition, the difference in the 
energy associated with the transition between the states '^l==n and ^'l=o is zero to leading 
order in N: for 7 —> 00, the occupation of some "quasi- momenta" states of quantum number 
m is a flat, Fermi-Dirac-like distribution, with \ifi\ < N/2 for L = [8]. To get to the 
state with L = iV -|- one has to promote the atom from the single-particle state with 

m = —N/2 to the one with m = A''/2 + 1. The difference between the energy of the two states 
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Fig. 3 - Schematic graphs showing AL in the occupancy of the quasi-momenta states with quantum 
number rh (high), and W in the dispersion relation £ = £{L) (low). 




is {N/2 + 1)^ — = + 1, which is indeed zero to leading order in N (which is iV^). 

For large but finite 7 the width of the Fermi-Dirac-like distribution within the Lieb-Liniger 
model - shown in the higher graph of Fig. 3 - is AL = N^/ (7 -f 2) + 0(1) [8], or 

N 

AL«iV-2— +0(1), (7) 

7 

neglecting terms of order N/j'^ [which is self-consistent with the final result, that implies 
that 7 ~ 0{N)]. If the single-particle operator AV acts on the state \'^l=n) once, similar 
arguments to the ones mentioned in the previous paragraph (for 7 — s- cx)) imply that the 
change of the total angular momentum which does not alter the energy of the system to 
leading order in N, is again equal to the width AL. 

However, in order for the transition to be allowed to first order in Ay, AL has to be 
larger than the width of the energy barrier W shown in the lower graph of Fig. 3 [defined as 
W = N — Lb, where £{L = Lj,) = £{L = TV)]. Since for large 7 the speed of sound for L = 
is Cs = d£/dL = iV7V(7 + 2f [9], therefore £(Lb) w £{<d) + c,Lh, and L^ = 1 + 0(1/7), or 

= AT- 1 + 0(1/7). (8) 

In order for W < AL (and the decay rate of persistent currents to be allowed to first order 
in AV), Eqs. lO and lO imply that the corresponding typical value of 7, denoted as 7decay, 
can be as large as N, [or gdccay can be as large as 0{N^)]. In other words, the decay rate 
of persistent currents is highly suppressed up to a value of .gdccay ~ 0{N'^), far beyond the 
coupling (7tg ~ 0{N) for which the system gets in the Tonks-Girardeau regime. Over the 
large window of couplings gxG ^ 9 ^ 5docay the system is very close to the Tonks-Girardeau 
limit and on the same time it supports persistent currents. 

To get numerical evidence about these estimates, we have calculated the pair correlation 
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function, defined as 

^ ^ ^ (A^-l)(*L=olE.'5(0-^.)l*L=o) ' ^' 

as function of the coupling g. Figure 2 shows 1 — 5^^''(0) (two higher curves in Fig. 2) for 
N = 3, and 6. The proximity of 1 — g^^\0) to zero/unity gives the extent to which the 
system is in the Gross-Pitaevskii/Tonks-Girardeau hmit. These results are consistent with 
the previous estimates gro ~ 0{N), and gdocay ~ 0{N'^). Furthermore, for a given g ~ 0{N), 
the matrix element decreases with increasing N, and probably approaches a step-like function 
for very large N, that reflects the Fermi-Dirac-like distribution. Figure 2 provides a specific 
example in this very small system, in which for 17x0 5; 5 5; Sdecay the system is close to the 
Tonks-Girardeau limit, while the decay rate of persistent currents is suppressed. 

Summary. - In conclusion, we considered bosonic atoms which rotate in a toroidal trap 
and examined the transition between the Gross-Pitaevskii and the Tonks-Girardeau limits. 
For any coupling £{L' ^ N - L) - £{L) ^ L' - L, and thus £{N) - £{0) = N. The energy 
barrier that the system has to overcome in order to get from the Z = 1 to the I = state 
and the corresponding distribution function of the atoms imply that persistent currents exist 
all the way between the Gross-Pitaevskii limit - above some critical coupling constant - up 
to a much higher critical coupling constant, close to the Tonks-Girardeau limit. Our model 
provides an example of a system where Bose-Einstein condensation is neither a necessary, nor 
a sufficient condition for superfluidity. 
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